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We extend the definition of the Jordan- Wigner transformation to three dimensions using the gen- 
eralization of ideas that were used in the two-dimensional case by one of the present authors. Under 
this transformation, the three-dimensional XY Hamiltonnian is transformed into a system of spinless 
fermions coupled to a gauge field with only two nonzero components. We calculate the flux per pla- 
quette for the 3 elementary perpendicular plaquettes of a cubic lattice, and find that it is nonzero for 
only two of the plaquettes. We provide a simple interpretation for the average phase-per-plaquette 
being 7r on the plaquettes where it is nonzero. Then we apply these new findings to the investigation 
of the Heisenberg bilayer antiferromagnet. 
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In 1928, Jordan and Wigner§ introduced a transfor- 
mation to write the spin operators S~ and S z in terms 
of Fermi operators in one dimension. For example, this 
transformation allows a mapping of the XY Hamiltonian 
in one dimension into a Hamiltonian of noninteracting 
spinless fermions. Trying to generalize this result to two 
dimensions proved to be difficult. But, a natural exten- 
sion of the Jordan- Wigner (JW) transformation was in- 
troduced in 1993 by one of the present authorsH to study 
the interchain-coupling effect on the one-dimensional 
spin-1/2 Heisenberg antiferromagnets. Several other au- 
thors have given implicit ways of defining a generalized 
JW transformation in two dimensions J3u The generaliza- 
tion of the JW transformation has been attempted even 
in three dimensions by Huerta and Zanellid who intra- 
duced an implicit transformation, and by KochmanskiEl 
whose transformation, however, does not preserve some 
of the spin-commutation relations. _ 

In this work, we follow the procedure of Ref.El to extend 
the definition of the JW transformation in three dimen- 
sions. This paper is organized as follows. We start by 
reviewing the JW transformation in one and two dimen- 
sions. Then we define its three dimensional version. For 
the XY model, we calculate the components of the ef- 
fective gauge field that results from the transformation, 
and the flux per plaquette as one spinless-fermion moves 
around it. Finally, we use this transformation to study 
the effect of the interlayer coupling on the Heisenberg 
bilayer. 

The JW transformation can be performed independent 
of any model Hamiltonian. In one dimension this is writ- 
ten as 
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with Q(i — j) = 1 for j < i and for j > i being the 
(discrete) Heaviside step function. Note that the phase 
term in this transformation is obtained by summing over 



all those sites to the left of site i. Also, when we calcu- 
late spin commutation relations we find that the fermion- 
number operator, at the site with the lowest index and 
hosting one spin of the commutator, is contained only 
once in the resulting phase (refer to Ref.el for more de- 
tails). This assures that the spin commutation relations 
are preserved. In two dimensions, a generalized transfor- 
mation based on this idea has been introduced in Ref . p] . 
This was written as: 
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where now a given site is specified by two indices, i along 
the x-axis and j along the y-axis. In this transforma- 
tion, the phase at site (i, j) is obtained by summing over 
all sites to the left of the vertical line passing through 
as well as all sites directly below This is a 

necessary and sufficient condition that allows the gener- 
alized JW transformation to preserve all spin commuta- 
tion relations. The phase term in @ can be written as: 

e^EHo^ 2 ^ 2 ' with 

w(x, z) = 0(xi - zi)(l - S XuZl ) 

+6(ar 2 - z 2 )S XltZl (l - 8 X2>Z2 ); (4) 



x = (xi, X2) and z = [zi, z 2 ) denoting sites on the square 
lattice. A necessary condition for the preservation of the 
spin commutation relation is that the following condition 
is satisfied 
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A relation that is indeed fulfilled by Eq. (^). 

Next , consider the transformation of the XY model 
under (|3|) , and calculate the flux per plaquette generated 
as one spinless fermion moves around an elementary pla- 
quette. The phase change is found to be: 



(6) 
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This result is more general than the one of Fradkinfl who 
reported that the phase should be 7T, and also shows that 
Shaofeng'sQ claim that the flux per plaquette is zero is 
not true. 

Using the same procedure described above for one and 
two dimensions, we define the JW transformation in three 
dimensions by: 
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where now the phase is: 
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x = (xi, X2, S3) and z = (zi,Z2, Z3) are now sites on the 
cubic lattice. It is straightforward to check again that 
Eq. (|^) is satisfied, and that all spin commutation rela- 
tions are preserved. To obtain the phase (||) we choose a 
plane going through the site where we want to define the 
transformation. Here the plane chosen is the one perpen- 
dicular to the vector (0,0, 1). Then we sum over those 
sites on this plane as in the two-dimensional transforma- 
tion. Finally we sum over all sites that are behind this 
plane. 

As for the phase change that occurs when a spinless 
fcrmion completes a motion around a plaquette, we con- 
sider the elementary plaquettes in the xy, xz, and yz 
planes. The flux per plaquette in the xy plane turns out 
to be the same as in two dimensions, that is 
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For a spinless fermion starting at site (i,j,k), rii t j t k = 1- 
Then for this fermion to go to site (i + l,j,k), this later 
site has to be empty due to Pauli exclusion principle; 
thus rii+i^k = 0. In this approximation one finds that 
the phase is n. For a plaquette in the yz plane the result 
is 



(11) 



Similarly, the phase per plaquette can be set to be ap- 
proximately 7r. As for the phase change around an el- 
ementary plaquette in the xz plane, it is found to be 
identically zero. 

The results for these fluxes per plaquette suggest that 
as a spinless fermion moves around a given plaquette, it 



will couple to a gauge field if the plaquette belongs to 
plane xy or yz, but no coupling to a gauge field occurs in 
the remaining plane xz. This thus suggests that the effec- 
tive magnetic field to which the spinless fermions couple 
has only two components. To show this, we consider the 
XY model written for coupled layers: 
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where J and J± are respectively the intralayer and inter- 
layer coupling constants. Under the JW transformation 
(0), the Hamiltonian ( |l2| ) gives 
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Here Vj^j.fc = </>i+i,j',fe _ 4>u 3 ,k and S7k<t>i,j,k = <f>ij,k+i - 
4>i.j,k designate the discrete derivatives along the x and 
z-axes, respectively. This effective Hamiltonian describes 
the motion of spinless fermions coupled to a gauge field. 
The later field is given by Ai = V 'i<f>ij t k along the x-axis, 
Aj = Vj4>i.j,k along the y-axis, and Ak = Wk<t>i,j,k along 
the z-axis. This leads to 
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Note that the y-component of A = (Ai, Aj, Ak) is found 
to be identically zero. The effective magnetic field is ob- 
tained through the discrete curl of A, that is B(i, j, k) — 
V x A(i,j.k), [vectors A and B depend on the position 
k)]. The components of B which are given by 



Bj = 0, 

B k = 7r(ni,j,fe - m+i t j t k), 



(15) 



are consistent with the results of the fluxes per plaquette 
in Eqs. ( |l(i| ) and (p"l|), (the unit area being 1). Therefore 
the effective magnetic field to which the spinless fermions 
are coupled has two nonzero components only. 

In the second part of this paper we apply the above 
findings to the study of the Heisenberg bilayer antiferro- 
magnet that has received a lot of attention in the past 
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few years.Blia The use of .the JW transformation and the 
bond-mean-field theoryBtJ is advantageous-, ie, compari- 
son with the Schwinger. boson approach ,HL3u3 modified 
spin-w avje a pproaches ,Ej or the bond-operator mean-field 
theory,! 1 tI 18 ! because no constraint on the number of par- 
ticle per site is required for the JW spinless fermions. We 
would likg_to draw the attention to the bond-mean field 
approachoO which has so far given very good results in 
the case of thti-Heisenberg ladder J2J and coupled Heisen- 
berg ladders.Eil The bilayer antiferromagnet is a system 
that presents an order-disorder second order quantum 
transition at zero temperature, and may be of some rel- 
evance to some of the high temperature superconductors 
in their low-doping normal state.E 2 ! 

The Heisenberg model on a bilayer reads as follows 

H= H X Y + J ^ S tjM S t+l,j,k 
i,j,k i,j,k 

Hamiltonian (|l6| ) is simplified using the approximation 
where the average flux per elementary plaquette is n on 
the xy and yz planes, and zero on the xz plane (as dis- 
cussed above). To achieve this, we choose the following 
configuration: the phases are alternated ...it — — tt — 0... 
along the adjacent bonds on the y-axis, and zero on all 
bonds on the remaining axes. Note that the alternated 
phases could be put rather on the x-axis without chang- 
ing the physical results because of gauge invariance. In 
the case of the bilayer system, the alternated phases can- 
not however be put on the z-axis because there is only 
one bond along that axis. Further simplification is done 
by using the bipartite character due to antiferromagnetic 
correlations. Then we decouple the interacting quartic 
terms by introducing the alternated magnetization pa- 
rameter mi = 2(Sl ] k ) = 2{c\,j,k c hj,k) -1 = (- l Y m , and 
the bond parameters Q = (ci,j,k4+i,j,k) = ( c *»i.fc c i,j+i,k) 
within the planes and P = (cij^cj . perpendicular 
to the planes. The Hamiltonian takes the following form 
in the reciprocal space: 

H=^2(2J+J±)m(^e lL -ftf k ) 

k 

+ ^[iJisinfc x +j(k y ,k z )}elfu + H.c. (17) 

k 

where ^(ky^kz) = Jicosky + J±iCosk z , = cj^ (resp. 
/k = c^) is the fermion operator on the sublattice A 
(resp. B). The quantities J\ and J±i are given by J\ = 
J (I + 2Q) and J±i = J±(l + 2P). The energy spectrum 
takes the expression: 

E±(k) = ±{(2J+ J^fm 2 

+J 2 sm 2 k x +j 2 (k y ,k z )} 1 / 2 (18) 

leading to the ground-state energy per site: 



E GS = ^-^m 2 + 2JQ 2 + J ± P 2 -J^E + (19) 

where by definition (d 3 k) = dk/(2ir) 3 . Because of the 
periodic boundary conditions along the z-axis, Ji is 
counted twice. Thus we have to divide the value of J± 
by 2 in all our equations, or multiply J± by two in the 
final results wherever Jj_ appears. 
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Interlayer coupling 



FIG. 1. The parameters m, Q and P are plotted as a func- 
tion of J±. Full lines are obtained for finite m, and dashed 
lines are for m = 0. The dashed-dotted line is E g defined by 
Eq. @. 

The parameters m, Q, and P are calculated self- 
consistently through the equations: 

m = J (d 3 k)m(2J + J±)/E+ 

Q^J {d 3 k ){ J^ 2 k x+l{ k yM co,k v} ,,E + 

P = J (d 3 k) 1 (k y ,k z )cosk z /2E+, (20) 

obtained by minimizing the ground-state energy with 
respect to m, Q and P. Two solutions for the mag- 
netization m/2 are possible depending on the value of 
J±, Fig. [[[ Let us first analyze the trivial solution 
m = of Eqs. (pif). In this case the energy spectrum 
E k = ijjj 2 sin 2 fcj, + 7 2 (fc y , fcz)} 1 / 2 presents an energy 
gap at k = (ko, n, 0) or (ko, 0, 7r), (ko — 0, it): 

E g - \J±i - Ji| = \J±{1 + 2P) - J(l + 2Q% (21) 

if J±i > J\. This condition is deduced from j(ky, k z ) = 
which is essential for a vanishing gap. Numerically we 
find that the value above which an energy gap opens 
is Jl = 2 x 0.78J = 1.56J, see Fig. |. However, 
this does not mean that this is the critical value for 
the order-disorder second-order quantum transition. As 
for the finite m solution, figure [j] shows that to ^ 
for Jl < J± c = 2 x 2.1J = 4.2J. The phase is or- 
dered antiferromagnetically for J± < J± c , whereas it is 
a disordered quantum state for J± > J± c . The value 
J± c = 4.2J compares wehjta.ihe result (4.48J) of the 
Schwinger boson approacli§E3lij-|and that (4.3) of the 
self-consistent spin-wave the£«y.E3 But it disagrees with 
the Monte Carlo simulationst3 result (2.5 J) or the series 
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expansionil calculation (2. 56 J). Chubokov and Morr0 
reported that while transverse spin-wave excitations are 
gapless (Goldstone modes) longitudinal spin-wave exci- 
tations are gapped for J± c > J± > 0. Our result seems 
to indicate that a finite value of the interlayer coupling 
(J_l = 1.56J) would be needed for longitudinal fluctua- 
tions to become relevant. Note however that we have to 
consider gaussian fluctuations about our mean-field point 
to account for the spin-wave excitations in the ordered 
phase. We expect that fluctuations beyond the mean- 
field point will lead to a better estimate of J± c . The mag- 
netization increases with Jj_, passes through a maximum, 
then decreases in agreement with previous results. til 

As for the parameters P and Q, when Jj_ increases 
P saturates at 0.5, whereas Q decreases to 0, Fig. [j]. 
For very large values of J± this approach becomes less 
accurate. Perturbations-expansion about the transverse 
dimers may be used.t3'E3 
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FIG. 2. The reduced ground-state energy is plotted as a 
function of the reduced interlayer cpupling as defined in the 
text. The dashed line is from Ref. fS\. 

Finally, the result of the calculation of the ground-state 
energy, Egs, is reported in Fig. || in order to support 
the validity of the present approach. To compare with 
the result of Weihongo obtained using Ising and dimcr 
expansions, v = Eqs/{^J + Jj_) is drawn as a function 
of 7/ = Jj_/(J + J_l). The maximum of v at rj = 0.5 
coincides with that reported by Weihong, and in general 
both results agree very well. 

In summary, we extended the definition of the JW 
transformation to three dimensions and gave a full anal- 
ysis of its consequences in the case of the XY model. 
This model is found to be transformed to a system of 
spinless fermions that are coupled to a gauge field with 
two components only (one of the three components be- 
ing zero). Accordingly, the flux per plaquette is nonzero 
for 2 elementary perpendicular plaquettes, but is zero 
for the third one. We gave a simple interpretation of 
the average flux being it for those plaquettes where it is 
nonzero. These findings were used to study the Heisen- 
berg bilayer antiferromagnet within the bond-mean field 
theory. We recalculated the order-disorder critical cou- 
pling and the ground-state energy. Very good agreement 
with previous results is found. Finally, the present three- 



dimensional JW transformation can be applied to models 
of three-dimensional quantum spin systems. For simplic- 
ity a cubic Bravais lattice was used in this work, but the 
transformation introduced here can be easily generalized 
to other Bravais lattices. 
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